
# Shepard-like diagram
# ********************

sphepardCophCor <- function(distMat, copheneticDist)
{
	plot(distMat, copheneticDist, xlab="Chord distance", 
	  #ylab="Cophenetic distance", asp=1, xlim=c(0,sqrt(2)), ylim=c(0,sqrt(2)), 
	  main=c("Linkage",paste("Cophenetic correlation ",
	  round(cor(distMat, copheneticDist),3))), )
	abline(0,1)
	lines(lowess(distMat, copheneticDist), col="red")
}
